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We explore the adiabatic particle excitations of an interacting field in a cosmological back- 
ground. By following the time-evolution of the quantum state corresponding to the particle 
excitation, we show how the basic properties characterizing the particle propagation can be 
recovered from the two-point propagators. As an application, we study the background- 
induced dissipative effects on the propagation of a two- level atom in an expanding universe. 



OO 

O I. INTRODUCTION 

o 

Much attention has been given in the literature to the different particle concepts existing 
in curved spacetime [1-3], which can be roughly classified into three classes. First, there 
are the global particle concepts which depend on special properties of the spacetime: if 
asymptotically flat regions exist, particles can be defined with respect to those asymptotic 
observers; if the spacetime has enough symmetries, particles can be defined with respect 
to them. Second, whenever there is a separation of scales between the propagating degrees 
of freedom and the background, particles can be defined in full generality, without making 
^jpj any reference to the specific form of the spacetime: these are the adiabatic and quasilocal 

particle concepts. Finally, the less restrictive particle concept is the operational definition 

,_i in terms of the response of a quantum mechanical detector. 

Most analysis have however dealt with non-interacting particles, i.e., with particles which 
are free except for the classical gravitational interaction with the curved background. When 
considering interacting processes in the universe, as for instance in nucleosyntehesis, the 

t-h curved background can be often safely neglected. However, there are situations in which 

\Q both the quantum interaction between the particles and the classical gravitational interaction 

with the background spacetime can be relevant: for instance it has been suggested [4-8] that 
this might be the case when considering the trans-Planckian problem [9-11] both in black 
hole physics and cosmology. 

While there are several works studying interacting particles in curved spacetimes (see for 

r> instance [1, 12-20] and references therein), most of them consider global particle concepts 

which rely on the existence of asymptotic regions or on specific symmetry properties of the 
spacetime. However, other particle interpretations which do not require particular properties 
of the spacetime can also be extended to the interacting case. In particular, when scales are 
well separated in cosmology, it is natural to look for an extension of the adiabatic particle 
concept. Indeed interacting adiabatic particles were already considered in Refs. [21, 22], 
which dealt with the disspative effects on the propagation a two-level atom in a an expanding 
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universe. 

In this paper we pursue the work of Refs. [21, 22] in two different directions. On the 
one hand we establish a general framework for the analysis of adiabatic interacting particles 
and quasiparticles in cosmology, thereby extending the results of Ref. [23] to cosmological 
backgrounds. We will construct and follow the time evolution of the quantum state corre- 
sponding to the (quasi)particle excitations, investigating how the two-point functions can 
be connected to observable quantities. The propagation of a quasiparticle in a physical 
medium and the propagation of a particle in a curved background have many similarities 
[24], and both situations can be treated using similar techniques. In particular, the closed 
time path (CTP) method [25-30] provides a generic way to deal with field theories over 
arbitrary spacetime backgrounds with an arbitrary field states. (For a brief introduction to 
the CTP method adapted to the notation of this paper see appendix A in Ref. [22].) 

On the other hand, using this framework, we concentrate on the novel dissipative effects 
induced by the universe expansion, beyond those present in flat spacetime. We expect 
those effects to be relevant when the interaction timescale is of the order of the expansion 
timescale. Indeed, in Ref. [22] it was mentioned that within this formalism one can readily 
show that the two-level atom becomes excited when propagating in a de Sitter spacetime in 
the presence of a conformal radiation field in the vacuum, the excitation rate corresponding 
to the effective de Sitter temperature [1] . In this paper we elaborate on this point by making 
explicit the calculation and generalizing the result to other backgrounds. 

It will be important to treat appropriately the several timescales appearing in the prob- 
lem: the propagation timescale E^ 1 , corresponding to the inverse de Broglie frequency of 
the particle; the interaction timescale t- mt ; the universe expansion timescale H -1 , given by 
the inverse Hubble rate, and finally the observation timescale i b s - Two timescale hierar- 
chies will be assumed thought the paper. First, in order for the adiabatic particle concept 
to be applicable, the typical propagation time E^ 1 must be much shorter than the typical 
expansion time of the universe H~ l , so that one can look for an analytic Wentzel-Kramers- 
Brillouin (WKB) approximation for the two-point propagators. Second, the observation 
timescale i b s will be assumed to be always much larger than the interaction timescale tint, 
so that the asymptotic field theory results can be (approximately) applied. Depending on 
the other timescale hierarchies several cases will be discussed in the paper. 

The paper is organized as follows. In Sect. II we present the relation between the propaga- 
tors and self-energies in curved spacetime, making the connection with the CTP formalism. 
In Sect. Ill we give a field-theoretic description of interacting adiabatic particles and quasi- 
particles in cosmological backgrounds by following the time-evolution of the expectation 
value of the Hamiltonian operator. In Sect. IV we apply the results to study the dissipative 
effects in the propagation of a two-level atom. Finally, in Sect. V we summarize and discuss 
the main results. 

Throughout the paper we use a signature (—,+,+,+) and a system of natural units with 
h = c = 1 . The same symbol is used for a quantity and its Fourier transform provided there 
is no danger of confusion. 
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II. INTERACTING FIELDS IN CURVED BACKGROUNDS 



In the following we consider the general situation in which there is a scalar field in a 
generic state p, over a globally hyperbolic spacetime characterized by some metric g^ u . Most 
work studying interacting fields over general curved backgrounds was developed during the 
late 70s and early 80s, and was focused on the study of the renormaliz ability of the theories 
(see Ref. [1] and references therein). More recent works have focused on general properties of 
interacting fields in curved backgrounds [31-35] and on the particular case of de Sitter and 
anti-de Sitter [36-40]. In this paper we will not attempt to make any review of the subject, 
nor make a complete presentation of interacting quantum fields in curved backgrounds; we 
will simply highlight some aspects relevant for us. 

The Feynman propagator, positive and negative Wightman functions and Dyson propa- 
gator, 



Gu(x,x') = Gf(x,x') 
Gi 2 {x,x') = G + (x,x') 
G2i(x,x') = G^(x,x') 
G 2 2(x,x') = Gb(x,x') 



Tr(pT4(x)$(x')), (la) 

■Tr(p4>{x)4>{x>)), (lb) 

■Tr(pj>(x'Mx)), (lc) 

Tr(pT0(x)0(z')), (Id) 



respectively, are correlation functions appearing in the CTP formalism (which is natural 
when dealing with interaction in curved spacetimes), and can be conveniently organized in 
a 2 x 2 matrix G a f,: 

n h(r ') _ ( G F (x,x') G-(x,x')\ 

u ab [x, x ) - ^ ^ j . w 

We may also consider the Pauli-Jordan or commutator propagator, 

G(x,x'):=Tt( P [^x)J(x')]), (3a) 

and the Hadamard or anticonmutator function 

GW(x,x') :=Tv( P {^x)J(x')}). (3b) 

Finally, one can also consider the retarded and advanced propagators, 

G R (x,x') := 6(x° - x'°)G(x,x') = 9(x° - x'°)Tt (p[4>(x),4>(x')}), (4a) 
G A (x,x') := 9(x'° - x°)G(x,x') = 9(x'° - x°)Tt (p[0(x),4>(x')}), (4b) 

Self-energies T, ab (x,x') are introduced through the Schwinger-Dyson equation: 
G ab (x,x') = G^ b \x,x') + J d^^^dVV=^Gic ) (^^)[-^ cd ^,/)]G' d6 (z',xO, 

(5) 

where G^ b \x, x') are the propagators of the corresponding free theory and an Einstein sum- 
mation convention has been assumed for the repeated CTP indices a,b,c. . . £ {1,2}. The 
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self-energy components T, ab (x, x') can be computed as the sum of all one-particle irreducible 
diagrams beginning in a vertex type a and ending in a vertex type b. 

It will be useful to work with rescaled fields, propagators and self-energies as follows: 

4>(x) := [-g{x)] l /\ (6a) 

G ab (x,x') := [- 5 (x)] 1 / 4 G a6 (x,x / )[-5(^)] 1/4 , (6b) 
£ ab (x,x') := [-g{x)\ l l^ ab (x,x')[-g{x')] l / A . (6c) 

With these definitions 4>(x), G ab (x,x') and T, ab (x,x') are (bi)scalar densities of weight 1/2. 
In terms of the bar quantities, the relation between propagators and self-energies becomes 
identical to flat spacetime: 

G ab (x,x') = G a ° b \x,x') + J d 4 zdV G^(x,z)[-it cd (z,z')]G db (z',x'). (7) 

The retarded propagator obeys a direct relation with the retarded self-energy Sr(x, x') := 
S n (x,x') + £ 12 (x,x'), 

G R (x,x') = G < £\x,x') + J d^z^z' G { l\x,z)[-it K {z,z')]G K {z' ,x'). (8) 

A similar relation holds between the advanced propagator and advanced self-energy. Another 
useful combination is the Hadamard self-energy, which is defined as T,^(x, x') = S 11 (x, x') + 
T? 2 {x,x') [or equivalently as S^^a;,^') = — T, 12 (x,x') — T, 21 (x,x')] and which is related to 
the Hadamard propagator through 1 

GW(x,x') = -iJ d 4 ydVG , R (x,y)s( 1 )(2/,y / )G A (y , ,x / ). (9) 

In a spatially flat Friedmann-Lemaitre-Robertson- Walker background, with metric 

ds 2 = -dt 2 + a 2 {t){dx 2 + dy 2 + dz 2 ) (10) 

the field operator can be expanded in conformal modes, <fi(t, x) = ^ k 4>k(t) e* k ' x . In practice 
this means that each conformal mode can be analyzed separately and that we can exploit 
conservation of the conformal momentum. Eq. (5) can be particularized to the cosmological 
case: 

G ab (t,t';k) = G ( f b \t,t';k)-i J ds ds' G<®(t, s; k)t cd (s, s'; k)G db (s', t'; k), (11) 

where the rescaled fields are 4>k(t) = a 3 ^ 2 (t)4>k(t). Two particular relations will be of interest: 
the one corresponding to the retarded propagator, 

G R (t,t';k) = G { °\t,t';k)-i f d S d S , Gf(t, S ;k)S R ( S , S ';k)G' R ( S , ,t , ;k), (12a) 



1 Eq. (9) is not valid in full generality, but requires remote initial conditions and a sufficiently dissipative 
behavior. See Ref. [41] for further details. 
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and the one corresponding to the Hadamard function: 

GW(M';k) = -i J dsds / G R (t,s;k)sW(s,s , ;k)G'A(s , ,t , ;k). (12b) 

In expanding universes the propagators, though space-translation invariant, are not time- 
translation invariant. Nevertheless, we can always express the propagator in a Fourier 
transform with respect to the difference variable A = t — t', while keeping T = {t + t')/2 
constant: 

G R (u, T; k) := J dA e^ A G R (T + A/2, T - A/2; k) . (13) 

In general, the Fourier transform does not help in simplifying the equations any further. 
However, if both typical propagation and interaction times are much shorter than the typical 
expansion time, i.e. if the hierarchy E^ 1 , t- m t <C t G bs *C i/ -1 . holds, the retarded propagator 
can be solved from Eq. (12a): 

G R (u,,T;k)« — = , (14) 

H4 0) (^,T;k)]-i+S R (cu,T;k) 

where the free propagator is approximated by [— iG^\u, T; k)] _1 ~ — lo 2 + k 2 /a 2 (T) + m? . 
For this expression to be valid for ultrarelativistic particles there is an additional condition 
on the momenta: |k|/a(T) <C {t 2 nt H)~ l . 

In principle, the perturbative evaluation of the propagators in curved spacetime can be 
performed in the same way as in Minkowski, taking into account the CTP doubling of the 
number of degrees of freedom. However, there are a couple of caveats. 

On the one hand, couple naive perturbation theory may be spoiled both because of ultra- 
violet and infrared divergences. While the leading ultraviolet divergences of the coupling pa- 
rameters of an interacting field theory in a curved spacetime are the same as the divergences 
of the same field theory in flat spacetime, there are subleading divergences which depend on 
the local curvature of the manifold. In general, in the spirit of effective field theories, one 
has to consider the most general Lagrangian which is compatible with the symmetries of 
the problem taking into account the gravitational background [42-44] . Infrared divergences 
may also arise. Calculations which are formally correct might be spoiled by the divergent 
infrared behavior of the theory, even if final results are not explicitly infrared-divergent. As 
opposed to ultraviolet divergences, which share a common structure for all theories, infrared 
divergences depend on the large scale structure of the spacetime, and therefore not many 
general things can be said about them. Infrared divergences usually appear when dealing 
with massless (or effectively massless) fields. Heuristically, one can make sure that the in- 
frared divergences do not play an important role by checking that the relevant contribution 
to all intermediate expressions is not governed by the far infrared modes. If this is the case, 
one can be reasonably confident of the infrared stability of the results; otherwise one has to 
deal with the infrared divergences in a case by case basis. 

On the other hand, in cosmological backgrounds even the computation of the free propa- 
gators is technically challenging. Exact closed analytic results can be obtained only in a few 
particular cases. The conformally coupled fields are one of this cases, in which the evalua- 
tion of the propagators reduces to the corresponding flat spacetime calculation through the 
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introduction of the conformal time coordinate. (We will encounter a conformally coupled 
field in Sect. IV.) Approximate closed analytic results can be also obtained under the adia- 
batic approximation. Whenever the energy of the mode is much larger than the expansion 
rate of the universe, the WKB solution provides a good representation for the propagator of 
the mode [1]. Thus, the adiabatic approximation, besides providing a well-defined particle 
interpretation, also facilitates the computation of the propagators. 

III. INTERACTING ADIABATIC PARTICLES IN COSMOLOGY 

In Ref. [23] the propagation of quasiparticles in physical media from a second quantized 
perspective was analyzed using both a real-time approach and a frequency-based approach. 
In this section we generalize the real-time analysis to include adiabatic particles propagating 
in cosmological backgrounds. In detail, we investigate how the basic features characterizing 
the (quasi)particle propagation in cosmology can be extracted from the two-point correlation 
functions. We will concentrate on the novel aspects introduced by the universe expansion, 
skipping most technical details, which can be found in Ref. [23]. 

In cosmology scalar particles 2 can be labeled by their conserved conformal momentum. 
When the energy of the particle is much higher than the expansion rate, adiabatic particles 
can be introduced. From a second quantized perspective, quasiparticles are characterized 
by their comoving momentum, their energy and their decay rate, this latter two quantities 
being time-dependent. At this point we should make clear that with "decay rate" we refer 
to the rate at which the probability of finding the particle at a given comoving momentum 
decreases. Therefore having a non- vanishing decay rate does not necessarily mean that 
the particle decays into a lower mass state: "decay" can simply means that the particle 
has changed its momentum. The motivation for this terminology is the second quantized 
perspective, wich implies that we focus on the mode corresponding to the particle momentum 
rather than on the particle itself. 

A. Quantum states and energy of the excitations 

For a start, let us consider the quantum states corresponding to particle excitations in a 
cosmological context. Since by assumption the observation times are much longer than the 
interaction times, the construction of the quantum states will be based on the asymptotic 
in-out representation of the interacting fields, conveniently adapted to the expansion of the 
universe. Notice that here "asymptotic" refers to the properties of the interacting propaga- 
tor, and not the cosmological model. Therefore we do not imply that the universe must have 
any asymptotic region: the asymptotic in-out representation of the fields is valid provided 
observation times are large enough, regardless of the details of the universe expansion. 

In the flat vacuum, the mode-decomposed interacting field operator, when acting on 
remote past or future times, and when evaluated inside a matrix element, can be approxi- 



2 For brevity, in the following we will not make the distinction between particles and quasiparticles, and will 
simply call "particle" every long-lived elementary excitation carrying momentum and energy, even if the 
background field is not necessarily in the adiabatic vacuum. 
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mated by a corresponding asymptotic free field operator [45, 46]. In turn, this free field can 
be decomposed in creation and annihilation operators of the asymptotic particle states of 
physical momentum p. Therefore the asymptotic representation for the field operator can 
be written <p p « [Z / {2E p )] l / 2 [a p + a_ p ], where E p is the physical energy of the particles 
and Z is proportional to the probability for the field operator to excite the vacuum with 
an energy E p . In Ref. [23] it is argued that an analogous construction can be extended to 
quasiparticle excitations in flat spacetime. Here we wish to point out that the representation 
of the field operator in terms of creation and annihilation operators can be straightforwardly 
extended to cosmology as a function of the conformal momentum k. The expression goes as 
follows: 

^ k «^^|y[a k (t) + aL k (t)], (is) 

where in this case E k (t) is the physical energy of the particle excitation at time t and 
Zk(t) is proportional 3 to the probability for the field operator to excite the background 
with an energy E^(t). It must be stressed that the above representation is an asymptotic 
relation, valid only when evaluated inside a matrix element in the large time limit. The 
parameter Z^(t) can be renormalized to one by rescaling the field; we will assume that such 
renormalization has been done in the following. 

In a flat spacetime, the one-particle state is created by the action of the aysmptotic field 
operator on the vacuum: |k) = ajjO). In Ref. [23] it was shown that when the vacuum 
is replaced by a slowly-varying background state p one can similarly find the one-particle 
state by acting with the creation operator on p. The further generalization to cosmology is 
straightforward: 

Pk'W ~ t- ^U(t,t )al(t )p(t )a k (t )U(t ,t), (16) 

7l k + 1 

where n k is the occupation number of the mode with momentum k (an adiabatic invariant), 
to is the initial creation time (assumed to be remote). The denominator ensures the proper 
normalization. The background state is assumed to evolve in timescales much longer than 
the interaction timescale. 

The time evolution of the expectation value of the Hamiltonian is given by 

£< + >(t,to;k) = TV [p£\t)H k (t)] =^rTr [p(t )a k (to)U(t , t)H k (t)U(t, t )4(*o)] ■ (17) 

where Hk(t) = i?k(£)[4(*)^k(^) + «Lk(i)a-k(i) + l] i s the Hamiltonian of the 2-mode ±k (see 
Refs. [23, 47] for a discussion of the convenience of including also the opposite momentum 
mode in general, although this is irrelevant in this presentation). Following similar steps 
as in Ref. [23], we apply the Wick theorem, assuming that the background state p(t) is 



3 The precise definition goes as follows. Z k (t)S(Ef 3 — E k (t)) = J2 a p Paa(t)\(a\<t> k \/3)\ 2 5(Ep — E a ), where |a) 
and \/3\) arc a complete set of orthornormal eigenvectors of the Hamiltonian spanning the Hilbert space. 
The density matrix p is assumed to be approximately diagonal in this basis. 
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approximately Gaussian, and find: 

£ (+) (Mo;k) « [p(t )a k (io)4(*o)] Tr [p(t )tf(*o,t)flk(fMMo)] 

+ £ k (i) Tr [p(t ) a k (t )t/(to, «)4(*M*> *o)] 
x Tr [p(t )C/(io, t)a*(t)U(t, t ) 4(*o)] }, 

Expressing the creation and annihilation operators in terms of the asymptotic field and its 
derivative, we arrive at the following suggestive expression: 



n k + 1 



2 v / E k (t)E k (t' 



G + (t,t';k) 



(18) 



where E^°\t) = Tr [p(to)£/(io 5 t)Hk(t)U(t, to)] is the energy of the unexcited mode. Thus, 
the energy of the mode corresponding to the particle as a function of time can be computed 
from the positive Wightman function. Recall that for the above equation to be valid the 
condition i b s S> tint has to be fulfilled, where i b s the observation time is t — t' . 

In the following we compute the evolution of the propagators [and hence the evolution of 
the expectation value of the energy, according to Eq. (18)] for several timescale hierarchies. 



B. Short observation times 

Under the assumption of observation times much shorter than the Hubble timescale, 
^obs *C H~ l , one can find a time representation of the propagator by Fourier-transforming 
back Eq. (14) . Since by hypothesis the observation times are also much larger than the 
interaction time, we can further approximate the self-energy by its value at the pole, similarly 
as in flat spacetime. Assuming a small decay rate, one finds [21, 41] 

G R (t, t'; k) = ^ sin [E k (T)(t - t')] e-r k CT)(*-t')/2 6{t _ ^ (19a) 

G«(M';k) = cos [EUT)(t-t')} e -|r k CT)(t-t')/2| (19b) 

with T= (t + t')/2, 

k 2 

El{t) :=m 2 + + Re£ R (£ k (t),t;k), (20a) 
a 2 (t) 

and 

r k (t) := --^Im£ R (E k (t),t;k). (20b) 

Introducing these expressions in Eq. (18) and taking into account that G + (t,t') = 
[G«(M') + Gr(M')]/ 2 for t > t', we obtain: 

E(t, H- k) = E^(t) + E k (t)(l + n k ) e -rm(*-*') . (2 1) 
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Therefore the energy of the quasiparticle excitation at time t is -Ek(t), and it decays in 
a timescale T(T). The factor 1 + n^, which is relevant only for non- vacuum states, is a 
consequence of the fact that slightly more than one quasiparticle is excited because of the 
Bose- Einstein statistics (see Ref. [23]). 

The short-time approximation is valid provided the following scale separation is verified: 
^k" 1 j tint ^ ^obs ^ H~ 1 . A careful analysis [21, 41] shows that for ultrarelativistic particles 
there is an additional condition on the momenta: |k|/a(T) <C (i 2 nt iT) _1 . 



C. Long observation times 



When the observation times are of the order of the expansion timescale we can no longer 
use Fourier-transform methods. Instead, we will look for a suitable WKB approximation 
for the propagators. Let us start by considering the equation of motion of the interacting 
propagator. By acting with the differential operator 



1 d 

aHt)dt 



(with R(t) being the Ricci scalar and £ being the conformal coupling parameter) on equation 
(12a), we get the equation of motion for the retarded propagator: 



1 d 

a 3 (t) di 



a3(t)|) + ^ + W) + -| y 1 



G R (t,t';k) 

a 3/2( t ) 3/2(f/) 



+ 



a 3 / 2 (t)a 3 / 2 (t>) 



dsE R (t,s;k)G R (s,t';k) 



(22) 



a3(t) 



S(t - 1') 



This equation so far is exact. Since we are interested in a first order adiabatic solution, we 
start by discarding all terms of the equation of motion which are of higher order in H/E^: 



Q t 2+ m + 2(j)j 



G R (t,t';k)+ f dst R (t,s;k)G R (s,t';k) = -iS(t-t') 



(23) 



The equation of motion for the retarded propagator in the rescaled time u = Ht is: 
d 2 o k 2 

U R {u,u;k) + — 

/no 



ou z a z {u) 



1 f u 

G R (u,u';k) + — / dv T lR (u,v;k)G R (v,u';k) = -iH6(u - u ) 

Expanding the non-local term as 

1 r u _ 
N(u,u') := — / dvT, R (u,v;k)G R (v,u';k) 

H Jun 



5El(u) + HT k (u) 

the equation of motion can be approximated by 

d 2 , x d 



du 



G R (u,u';k) + 0(H 2 ), 



(24) 



H2 d^ 2+HT ^d-u +E ^ u) 



G R (u, u; k) = —iH5(u — u ), 



(25) 
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where E^(u) = m 2 + k 2 /a 2 (u) + 5E^(u). The leading order adiabatic solution is (assuming 
that I\ is much smaller than E^) 

G R (u,u';k) = - ~* =sinf ^- P dv E k (v)] e~ *n D d " r k(") 0( u - u), 
or, going back to the original time t, 

eK(M ' k) - vrnmrr (/,' ds£tW ) - o- <*> 

The interacting propagator can be expressed in terms of the adiabatic evolution of the 
quasiparticle energy and decay rates. Similarly, the Hadamard propagator can be expressed 
as [41]: 

e,1,(f '^ k ^7&^^(/' ds H eHJ;id ' rt, ' ,/21 - (27) 

Taking into account that G+(t,t') = [G R (t,t') + t')}/2 for t > t' , and introducing 

Eqs. (26) and (27) in (18) we find in the leading adiabatic order: 

£ (+) (i, t ; k) « £7(°>(t) + E k (t)(l + n k ) e~ # dsF « . (28) 

From this equation we identify E^t) as the energy of the quasiparticles, which is slowly 
evolving with the background, and T^t) as their net decay rate. Recall that the factor 
(1 + n k) is a consequence of the fact that the quasiparticle state actually contains more than 
one particle excitation. Recall also that in this context decay does not necessarily mean the 
particle going into a lower mass state, but it means that the particle changes its momentum. 

Let us see now how the energy and the decay rates can be extracted from the propagators, 
by considering separately two different situations. 

a. Short interaction times: E^ 1 ,^ <C t bs,^ _1 - When the interaction times are 
much shorter than the inverse expansion rate the particle energy and the decay rate can be 
computed from Eqs. (20a) and (20b), respectively [21, 41]. In other words, as far as the 
interaction process is considered the universe expansion can be neglected. 

For non-relativistic particles the requirement that the interaction time is small as com- 
pared to the typical expansion time can be somewhat relaxed, since the precise condition 
for Eqs. (20a) and (20b) to be valid is: 

t mt dE k k 2 /a 2 . . 

b. Long interaction times: E^ 1 <C tint,// -1 <C t bs- When the interaction times are of 
the order of the expansion timescale H -1 , the expansion effects cannot be neglected during 
the interaction time, and the Fourier transform is no longer useful. To investigate how SE^ 
and Tk can be expressed in terms of the self-energy, we reconsider the non-local term: 

N(u,u') = - T f U dv 7 SR( "^ ; E sinfl f V dv'E k (v')) e~G dv ' 9{v _ 
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Splitting the interval [u',v] in [u,v] minus [u,u'] yields 

N(u, u) = - f H„ -^ U ^^ L P S:^v'T^v')/(2H) 

sin J dv' E k (v')\ cos J dv E^{v') 

^ J\v' E^v'^sin(^ J\v' E^v'^j . 



— cos 



Replacing v by u in the argument of the exponential function (the total amount of decay 
during the interaction process is negligible; otherwise the particle would have completely 
decayed at the observation point), we find: 

N(u, u') » (5E k )\u)G R (u, k) + HT k (n) 9G ^< k ) , (30) 
where in non-rescaled time t, 



(31b) 



Eqs. (31) give the most general representation of the energy shift and the decay rate in 
terms of the self-energy. Notice that the self-energy is evaluated in a kind of frequency 
representation, but with the frequency varying along the interaction range, because the on- 
shell position changes significantly during the interaction process. Notice also the presence 
of the square root prefactor, which can be interpreted as the geometric mean of the particle 
redshift along the interaction time. 

Let us concentrate on the dissipative effects in the remaining of this section. When 
the interaction times are comparable to the inverse expansion rate, the interpretation of 
Eq. (31b) as a decay rate may be hindered by the universe expansion. In effect, since that 
the decay rate is only a meaningful concept when considering observation times much larger 
than the interaction time, the decay rate is also expected to evolve on scales much larger 
than the typical expansion time. However, the decay rate extracted from Eq. (31b) evolves 
on timescales comparable to the universe expansion because of the redshifting. To avoid 
these difficulties, let us consider the decay rate in the particle rest frame, which naturally 
can be conjectured to be: 

7k(t) = ~J* ds s ; k ) sin ds ' ^k(s')) • (32) 

Recalling that this quantity is only meaningful when observed for large periods of time, let 
us average the decay rate for timescales much smaller than the observation time. In other 
words, we choose a timescale t av much larger than ti nt but still much smaller than t ^ s , and 
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average the decay rate over that scale: 

rt+tw/'Z rv I rt" 

/ dt" ds£ R (i",s;k)sin / ds' E k (s' 

Jt-ta.v/2 Jt' \Js 



7k(t)«— / dt" 7k (i") 



| /-t+tav/2 ,■■/" / /•*" \ ( 33 ) 



"T-^av Jt-tw/2 Jt' 

Taking into account that t av is much larger than the interaction time, and introducing the 
semisum and semidifference coordinates we get: 

i /-t+tav/2 r ( rT+A/2 \ 

7k(t) ~ — / dT / dA £ R (T + A/2, T - A/2; k) sin / ds' E k (s') , 

mt av Jt-tav/2 J \JT- A/2 / 

or, taking into account that Er does not depend significantly on T during the averaging 
time, 

7k(t) ~ Im / dA S R (t + A/2, t - A/2; k) e 1 J *-a/2 ds ) . (34) 

777 J 

In terms of the following "improved" frequency representation for the self-energy, 

£ R ([£ k ],t;k) := y dAS R (t + A/2,f-A/2;k)e ! C 2ds£k W (35) 

the decay rate in the particle rest frame can be expressed as 

7k(i) ~ -^ImS R ([E k ],t;k). (36) 

Notice that the improved frequency representation has been defined only for on-shell values 
of the frequency; in this sense it cannot be considered a proper integral transform. Note 
also that the improved frequency representation yields the same results as the standard 
frequency representation for comoving particles. 



IV. EXPANSION-INDUCED DECAY IN COSMOLOGY 



Let us apply the techniques shown in the previous section to analyze in a particular 
model the dissipative effects on the propagation in an expanding universe. Let us consider 
two massive fields with large masses but with a small mass difference Am := M — m <C M. 
As it is shown in Refs. [48, 49], the model can be interpreted as a field-theory description of 
a relativistic two-level atom (of rest mass m and energy gap Am) interacting with a scalar 
massless radiation field. The masses of the fields will be assumed to be much larger than the 
expansion rate of the universe, allowing us to introduce the adiabatic approximation. The 
mass gap between the two massive states will be taken to be of the order of the expansion 
rate of the universe. 
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k-q 



FIG. 1: Feynman diagram leading to the one- loop self-energy of the less massive particle. Straight 
lines represent the fundamental sate, double lines the excited state and curly lines the masslcss 
particle. 

In detail, the model consists of two massive fields 4> m , and 4>m, interacting with a massless 
field, X) y i a a trilinear coupling. The total action can be decomposed as S = S m + Sm + 
S x + Si n t, where: 

S m = \J dtd 3 xa 3 (t) {{dt^f - J^(dM 2 - m 2 <£^ , (37a) 
Sm = \J dtd 3 xa 3 (t) {{dt^M? ~ ^{d^M? ~ M 2 4> 2 M ^j , (37b) 
S x = \f dt d 3 xa 3 (i) {{dtx? - J^(dxX) 2 ~ W)x 2 ) , (37c) 
Sint =gM j dtd 3 xa 3 {t)4> m 4>MX, (37d) 

with R(t) being the Ricci scalar. We assume that the massless field is conformally coupled 
to gravity, so that £ = 1/6. In this paper we will be interested in studying the vacuum 
effects, so that the massless field will be assumed to be in the conformal vacuum and the 
massive fields in the adiabatic vacuum.; see Refs. [21, 49] for an analysis of the temperature 
effects. 

The goal is to extract the dissipative effects on the propagation of the unexcited atom from 
the self-energy. To one loop, the self-energy can be computed by using CTP perturbation 
theory as (see figure 1) 

-iE^tajk) = {tgM) 2 c aa c bb j ^&<°\t u t 2 ; k - q^ft, t 2 ; q) (38) 

(no summation implied), where G*J®\ti,t2; k — q) are the free adiabatic propagators of the 
more massive field and A ( ^(t 1 ,t 2 ; q) are the free conformal propagators of the massless 
field [1, 21]. We recall that a,b,c... indices refer to the different CTP branches and that 
c ab = diag(l,— 1). The retarded self-energy corresponds to XR,(ti,t2;k) = S 11 (t 1 ,t 2 ;k) — 
S 12 (ii, t 2 ; k). Since we are investigating the effects induced by the universe expansion, and 
therefore the interaction time is of the order of the inverse Hubble rate, we will evaluate the 
self-energy in the improved frequency representation. 
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FIG. 2: Feynman diagram contributing to the decay rate of the field <p m . The diagram violates 
energy conservation, and it is hence forbidden in flat spacctimc. However this needs not be the case 
in an expanding universe. 



The imaginary part of the self-energy is given by 

lmt R (h,t 2 ) = 1 [X 12 (ti,t 2 ) - Z 21 (h,t 2 )] . (39) 

By applying the CTP Feynman rules in the physical time representation to the diagram in 
figure 1, to order g 2 we get: 

ig 2 M 2 f d 3 q e^g dt ' E t-^ti d *'N 1 
S {hMk) ~ a( tl )a(t 2 ) J (2vr)3 2 (t2 ) 2|q| ^ 



with (t) = (k-q) 2 /a 2 (t)+M 2 . Eq. (40) can be interpreted in terms of the square of the 
Feynman diagram represented in figure 2, corresponding to the emission of a photon (this 
diagram would vanish on shell in flat spacetime). In the (standard or improved) frequency 
representations, 

£ 21 ([£ k ],T;k) « -2iIm£ R ([£ k ],T;k). (41) 

The equality is valid because the Fourier transform of S 12 gives exponentially damped 
contributions in the high mass limit we are using. Therefore 

r*i a+i i 



ImS R( [^ k ],T;k) = - 5 2 M 2 / - d f / d * J- 

7 2o(ti)a(t 2 )7 (2vr) 3 2 A / J B*_ q (t 1 ) J B*_ q (t 2 ) 



(42) 

where ii = T+A/2 and t 2 = T — A/2. Noticing that it vanishes in the case of flat spacetime, 
let us evaluate this quantity for two different models of the universe expansion. 



A. Vacuum effects in de Sitter 

Let us consider flat coordinates in the de Sitter spacetime, which corresponds to a flat 
cosmological model with 

a(t) = o(T) e H{t - T ^ . (43) 



4 Notice that here we are abusing the notation because in the time representation the retarded self-energy 
is purely real. "Imaginary part" here refers to the standard or improved frequency representations. In the 
time representation, it corresponds to the odd part of Sr under the exchange of t\ and ti. 
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For the sake of simplicity, we will show the details of the calculation in the comoving case, 
and simply quote the result for the case of arbitrary geodesic motion. 

From Eq. (42), the vacuum contribution to the imaginary part of the self-energy in de 
Sitter is given by: 



//>oo 
dA J dk 



g 2 mk 



x exp 



167r 2 a 2 (T) 

T+A/2 



-i dt ( Am + -4-r e 11 ^ 



(44) 



/T-A/2 

Integrating the argument of the exponential yields, 

t v < t 1 n\ (a\ f°° ai 9 2 mk ^ ro A--Ju sin h(4£)-efc 

ImS R (m,T;0) = -y dAy dk 16 ^ 2q2(t) e ^ j 

where have added a small — ek term to ensure convergence of the k integral: 

g 2 mH 2 e (- iAm +V A 



r a 2 mH 2 e (-iAm+n)A 



2 ' 



The A integration can calculated by residues in the complex plane, closing the circuit through 
a parallel line passing by 2-rri/H and taking into account the pole at A = it. The result is 
the following: 

— q 2 vnAvn g 2 

ImS R (m,T;0) = ~^ e2wAm/H ^ = -—mAmn H/27r (Am), (45) 

where ut{x) is the Bose-Einstein function corresponding to a physical temperature T. This 
results coincides with the imaginary part of the self-energy in a Minkowski thermal bath 
at physical temperature H/(2tt) [21, 49]. A decay rate can be associated to this imaginary 
part of the self-energy: 

1 - Q 2 
T (T) = 7 o(T) = Im£ R (m,T;0) = ^Amn H/27T (Am). (46) 

Recall that in the comoving case there is no difference between the standard and improved 
self-energies. 

The result is not at all unexpected, since it is well known [1] that a comoving particle 
detector in de Sitter conformal vacuum perceives a bath of thermal radiation at a temper- 
ature H/(2tt). The fact that the decay rate has an exact thermal character can be traced 
to the existence of thermalization theorems under the presence of event horizons [50] . Our 
results are also consistent with those of Bros et al. [20] , who similarly find inest abilities for 
massive particles propagating in de Sitter, using a global rather than an adiabatic approach 
to the particle concept. 

This result can be extended to moving particles with arbitrary momentum. Notice that in 
this case it is important to use the improved frequency representation because the interaction 
timescale is of order H~ l . We have evaluated the integrals in a non-relativistic expansion in 
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k 2 /m 2 . Interestingly, we have found that the additional momentum-dependent terms cancel 
and that the result is 

ImS R ([£ k j,T;k) = ^-mAmn H/27T (Am) (47) 

at least to order k 10 /m 10 . It is therefore reasonable to assume that the above result is valid 
to all orders in perturbation theory. The result is not unexpected, since the decay rate in 
the particle rest frame 

1 ~ q 2 
7 = i m £ R ([£ k ] , T; k) = n H/2lT (Am) (48) 

is time- and momentum- independent as required by the de Sitter invariance of the problem 
(bear in mind that both the adiabatic and conformal vacua are de Sitter invariant). Since 
the invariance is not manifest from our expressions, this provides a non-trivial check of the 
calculation. We have also computed the decay rate in the comoving frame, according to 
Eq. (31). The result is the following: 

r k (t) = ^Amn H/2n (A m ) (l - ^^2) +0(k 4 /m 4 ). (49) 



B. Vacuum effects in power-law inflation 



Let us now consider the case of power-law inflation, where the scale factor evolves ac- 
cording to 



a(t) = a(T) 



(50) 



The Hubble rate is given by H(t) = a{t)/a{t) = a/t. Since expressions rapidly become 
cumbersome, we shall only sketch the calculation in the simplest situations, and will present 
the result for the other cases. 

As before, let us first consider comoving particles. According to Eq. (42), the imaginary 
part of the self-energy is given by 



ImS R (m,T;0) = - J dA J 



/0 ° „ 9 2 mk 
ak 



167r 2 a 2 (T) 

T+A/2 



1 



x exp ^ — 1 l dt 

'T-A/2 



Am + 



4T 2 
k (T 



o(T) V t 



(51) 



Integrating the argument of the exponential, and performing the integral over k yields (it is 
necessary to add a — ek term to ensure convergence) 



ImS R (m,T;0) = 



2g 2 m(-l + a) 



/dA 



-iAmA 



(2T + A) 



2a 



(52) 



{T a (2T + A) - (2T - A) 1 - Q [T(2T + A)] Q } 2 ' 
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In order to proceed further a particular value of a has to be chosen. Closed analytic expres- 
sions can be obtained for positive integers values of a. For the case of a = 4 we obtain: 

9g 2 m f e~ iAmA (-4T 2 + A 2 ) 4 

Im S R m,T; 0) = -f-^ / dA ^ 9 ' , a = 4. 

V 7 16vr 2 7 A 2 (2T 2 + A 2 ) 2 

This function can be integrated in the complex plane by considering the residues of the 
double pole located at A = —2-/3iT. 5 The result is the following: 

ImS R (m,T;0) = -^mAm e" 2 ^ mT = -^mAm e~ 8 ^ Am / H ( T ), a = 4. (53) 

This expression can be generalized, on the one hand, by considering different values of 
a, and on the other hand by considering particles in motion. We have found that at this 
order there is only a non-vanishing contribution to the imaginary part of the self-energy for 
even values of a larger than two. In the case a = 6 the following result is found: 

Im S R (m, T; 0) = — (3 + ^ 9 - mAm e -^V^UEAmT ? Q = g (M) 

Expressions become increasingly cumbersome with a. The momentum corrections can be 
found in Ref. [41]. 

The decay rates can be immediately extracted from the imaginary part of the self-energy 
following Eq. (36). Qualitatively, the behavior is the same as in de Sitter: the vacuum- 
induced decay rate is exponentially suppressed except when H > Am. However, in this case 
the rates do not have an exact thermal character and depend on time, because the Hubble 
rate varies during the universe evolution. Unfortunately we have not been able to compute 
the decay rate in the ultrarelativistic limit, which is the most interesting limit (because even 
slowly moving particles will become rapidly moving with respect the comoving frame as we 
go back in time along the universe expansion). 



V. SUMMARY AND DISCUSSION 

In this paper we have analyzed the adiabatic particle excitations of an interacting field 
in a cosmological background. In particular, we have shown that the time evolution of 
the energy of interacting adiabatic particles can be expressed in terms of the two-point 
correlation functions [Eq. (18)]. These in turn can be calculated in a WKB approximation — 
see Eqs. (26). The long-time dynamics of the particle can be essentially characterized from 
a second quantized perspective by two semi-local quantities, the energy shift and the net 
decay rate, this latter quantity being understood as the total transition rate to any other 
field state. 

Two different physical situations have been distinguished. If the interaction timescale is 
much smaller than the typical expansion rate, the energy shift and the decay rate can be 



5 The other double pole would be located at A = ie had we kept track of the e terms, and therefore does 
not have to be taken into account. 
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extracted from the retarded self-energy in a quasilocal frequency representation [Eqs. (20)]. 
If, on the contrary, the interaction timescale is comparable to the expansion rate, then 
the energy shift and the decay rate can be extracted from a kind of frequency transform 
which takes into account the evolution of the on-shell condition during the interaction time 
[Eq. (31)]. We also defined an improved frequency representation, Eq. (35), from which the 
decay rate in the particle rest frame can be extracted. 

These results have been applied to the analysis of the evolution of a doublet of massive 
fields interacting with a third massless field. We have focused on the dissipative effects 
generated by the universe expansion, which are relevant when the expansion rate of the 
universe is of the order of the mass gap or larger. In the case of de Sitter, the decay rate 
corresponds to that of a particle at rest in a flat thermal bath at the de Sitter effective 
temperature [Eq. (48)]; moreover we verified that the decay rate is de Sitter invariant. In 
the case of power-law inflation, a qualitatively similar decay rate is found [Eq. (53)]. 

The approach presented in this paper is complementary to other analysis of the cosmo- 
logical decay rates in the literature [13, 19, 20, 51], which have been based on asymptotic or 
symmetry-based particle concepts and S-matrix computations. Concerning the operational 
definition of the particle concept in terms of the response of a particle detector, in Ref. [41] 
it is argued that the operational and adiabatic approaches lead to identical results within 
the range of validity of the adiabatic approximation, even if self-interaction is taken into 
account. 

A key point in this paper has been using several approximations controlled by differ- 
ent expansion parameters. The adiabatic approximation, controlled by the ratio H/E^, 
provided us with a well defined particle concept and furthermore allowed having simple an- 
alytic expressions. The asymptotic approximation, controlled by the ratio iintAobsj allows 
to obtain model-independent predictions on the evolution of the propagators. The Gaussian 
approximation, which is formally controlled by 1/N (with N being the number of fields), 
reduces the problem of computing the evolution of the quantum state of the particle to the 
computation of the two-point correlation functions. Most of the times we have not explicitly 
indicated the expansion parameter of the different approximations for the sake of brevity. 

One of the requirements of the calculation of the dissipative effects in the adiabatic vac- 
uum is that observation times must be much larger than interaction times. Since interaction 
times are of the order of the Hubble time H~ l , interesting results can only be obtained when 
the Hubble rate is slowly varying, or, in other words, when H <C H 2 . This condition can 
only be achieved in inflationary contexts. Notice that in the power-law inflationary model 
H = H 2 /a, so that large values of a should be considered (although for computational sim- 
plicity we just studied the a = 4,6 cases). Concerning also this inflationary model, recall 
that no contribution to the decay rate was found for odd values of a; this can be probably 
attributed to some symmetry of this particular case, and most likely does not extend to 
higher order calculations. 

Let us end by mentioning that the decay rate, derived from the imaginary part of the 
self-energy, has a secular character, as expected. Even small decay rates could give an im- 
portant effect when integrated over large periods of time. Moreover, dissipation is a generic 
phenomenon which appears even in the vacuum. Therefore, if one considers remote times 
or large time lapses, the dissipative processes will become relevant, specially in situations in 
which the expansion rate, or the curvature, are large. 
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